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ABSTRACT 

A  new  and  useful  parameter  (a)  for  avalanche  diodes  is  obtained 
which  possesses  the  properties  associated  with  negative  conduc¬ 
tivity.  It  is  shown  how  <7  unifies  the  description  of  various  as¬ 
pects  of  device  behavior  such  as  diode  impedance  Z,  total  current, 
and  the  effect  of  device  radius  on  performance.  A  greatly  simpli¬ 
fied,  approximate  formula  for  Z  is  obtained,  in  terms  of  a,  which 
predicts  reasonably  well  the  significant  trends,  zero-crossings, 
and  peaks. 


Accepted  for  the  Air  Force 

1  ranklin  C  .  Hudson 
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NEGATIVE  CONDUCTIVITY  IN  SOLID  STATE 
AVALANCHE  D  1  O  D  K S 


i.  1  NT  ROD  I JCTION 

The  resistance  R  of  a  simple,  classical  resistor  is  calculated  in  an  elementary  fashion 
from 


which  relates  the  geometric  factors  (L  length,  A  cross-sectional  area)  and  maternal  param¬ 
eter  (o*  conductivity)  to  R.  The  impedance  Z  of  a  simple  parallel-plate  capacitor*  of  separa¬ 
tion  I  ,  and  cross-sectional  area  A,  completely  filled  with  a  lossy  dielectric  of  permittivity  < 

and  conductivity  a  ,  is 
o 

(jwf  +  a  )  A  *  u 

o 

This  r  eport  considers  the  description  of  an  avalanche  diode  in  a  new  simplified  manner  using  an 
effective  conductivity  on  The  impedance  Za  of  the  avalanche  region  will  be  shown  to  be  given 
approximately  by 

y  -  L  ,  , 

a  (jwt  4^)  A  * 

In  addition,  the  total  current  density  (particle  plus  displacement)  will  be  shown  to  be  specified 
without  approximation  by 

•Lp  (cr  +  juu  )  B  ^  (  D 

(where  B^  is  an  amplitude  constant)  and  the  radial  variation  of  the  elect  r*ic  field  and  currents 
will  be  shown  to  be  given  for  the  normal  IMPATT  mode  by  J  (Tr)f  where 

2 

T  -Jo;[io(ct1  +  )  (^) 

ui  material  per*meabi lit v ,  and  .!  (x)  is  the  Bessel  function  of  the  1st  kind  and  of  zeroth  order 

Previous  "exact”  expressions  for*  Z  in  the  literature  at*e  algebraically  complex  expressions 

which  require  lengthy  numerical  calculations  before  their*  content  can  be  made1  explicit.  The 

dependence  of  the  total  current,  1^,,  on  radian  frequency  (w),  IK*  current  density  bias  (  !(}).  and 

the  material  parameters  has  not  been  previously  developed  in  the  literature.  Finally,  the  radial 

variation  of  1  I  qV  n,  and  .1  — q\  p  (the  electron  and  hole  current  densities  with  \ 

x  n  M  s  p  s'  s 

saturated  cl r* i ft  speed)  has  not  been  previously  evaluated,  although  they  may  begin  to  be  noticeable 
in  sufficiently  large-diameter  ring  diodes  of  the  type  described  by  Gibbons  anil  Ylisawa.* 

11.  TOTAL  (’ LR RENT  DENSITY 

The  total  current  density  is  defined  by 

•  1 -  .1  t-  .!  +  ju?  c  E  .  ( ('i ) 

I  n  p  *’  x 

It  is  shown  in  the  Appendix  that 


1 


(7) 


e  =  y  e. 

X  Li  1 

i  1 


(  'jKix\ 

(E^e  *)  , 


j  +  j  y  a.E. 

n  p  Li  ii 


(8) 


i=  1 


where  B^,  B^  are  amplitude  coefficients,  and 


CTi 


—  Z ic oa' J  V 
J  o  o  s 


o>2  -  (K.V  )2  +  2jc oa  V 
is  J  os 


(i  =  1,  2,  3) 


(9) 


are  temporally  and  spatially  dependent  conductivities,  =  0,  and 


K0  =  -K,  /(w/V  )2  +  2j(w/V  )  a  -2(vM  /e  Vo 

Z  3  v  '  s  J  '  s  o  o  o'  s 

are  wavenumbers.  Here  a  and  n  1  are  the  ionization  coefficient  and  dcv/dE  evaluated  at  the  DC 

o  o 

electric  field  value  Eq.  By  direct  addition  of  (8)  and  jcoe  times  (7), 


Jrji  =  (a,j  +jwe)B^  +  ZB^cr^  +  jcoe  )  cos  (K^x) 
Direct  calculation  yields  =  —  jc^e,  which  leaves  the  final  result 

J  rp  ((7  j  "t"  )  8  j  * 


where 


Zio’J  V 
J  o  o  s 

u>  +  Zifv  V 
J  o  s 


:  Ws\  B 

V  ^  +  2jQ'oVs  /  1 


(10) 


(11) 


(12) 


for  which  normalized  values  are  plotted  in 
Fig.  1. 


Fig.  1.  The  real  and  imaginary  components 
of  the  total  current  density  (.Bp)  as  a  function 
of  frequency  (f )  for  an  X -band  silicon  diode  in 
which  the  avalanche  zone  is  1.5  microns  in 
length,  .Tn  lO^amp'cm^,  q,  6.6  x  10^, 
0.164,  and  «  lO-*2. 
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III.  D  E  V I C  K  1 M  P  E  D  A  N  C  E 


The  impedance  is  defined  by 

Za  ,  V/JtA  ,  (13) 

where  A  c  ross-sectional  area,  and 

r  +  L/2 

V  *  -  \  K  dx  (14) 

J-L/2  x 

is  the  voltage  across  the  diode.  Neglecting  space-charge,  integration  of  Eq.  (6)  (see  Appendix) 
gives 

E  constant  B. 

x  1 

so  that 

V  -  -E  L  -  —  B ,  L  ,  (IS) 

x  1 

I  rom  Eqs.(ll),  (13),  and  (15),  Eq,(3)  is  obtained  for  Z  .  The  exact  result  (see  Appendix)  is 


Z 

a 


/  IV)  s  in  0?  \ 
) 


( 1M 


where 


0,  _  K-L/2  (L/2)  /(w/V  )2  +  2j(w/\  )  <v  •  2«\1  /(V 

2  2  \f  's  s  o  o  o  s 

and 


~  <  (u>  t  2jn  \  )  (jo;  cos0?  -  V  K-  sin07) 

Figures  2  through  5  show  that  the  activity  threshold  and  peak  impedance  frequencies  air 
reasonably  well  predicted  by  the  approximate  result  of  Eq.  (3)  (in  conjunction  with  the  drift  /.one 
impedance),  although  the  values  of  in  the  vicinity  of  these  points  are  not  precise.  Thus,  the 
significant  trends,  peaks,  and  zero-crossings  are  revealed  by  Eq.(3),  although  the  magnitudes 
show  only  semi-quantitative  agreement  in  certain  ranges.  However,  this  is  a  far  better  approx¬ 
imation  for  than  provided  by  previous  simplified  analyses,  such  as  that  of  Gilden  and  Ilines,^ 
which,  for  example,  shows  no  negative  resistance  effects  associated  with  the  avalanche  region. 


IV.  RADIAL  VARIATIONS 

A  lengthier  analysis^  reveals  that  the  usual  quasi-static,  one-dimensional  picture  may.  in 
a  more  accurate  field  theory,  be  described  as  a  quasi -T  EM  radial  wave  mode.  One  may  an  (’nipt 
to  approximate  this  by  a  TEM  radial  wave  mode  (which  wdll  only  satisfy  boundary  conditions 
approximately).  The  equations  for  a  radial  mode,  with  the  only  non-zero  components  being 
and  I!  ,  are 

ip 


3 


(18) 


i  a 


-  (rH  .)  =  =  (d.  +  jw€  )  E. 

r  9r  (pi  Ti  l  J  1 


ar 


ium  H  . 
J  o  01 


(19) 


which  combine  to  yield 


a  e 

—  rf-g— =  jo;p  (cr.  +  joK )  E.  =  T.2E. 
r  ar  \  9r  /  J  o  i  J  i  l  i 


The  non-singular  solution  to  Eq.  (20)  is 


Ei  =  Jo<V>  Bi  e 


-jK.x 


(i  =  1,  2,  3) 


(20) 


(21) 


By  direct  calculation  T,  — jco  |jl  (cr.  +  j  cot),  shown  in  Fig.  6,  while  a  -  <7_  -jun:  implies  that 

1  O  1  Lb 

T  T  0.  This  behavior  is  peculiar  in  that  the  radial  TM  wave  mode,  with  non-zero  E  ,, 

2  3  x 

H  and  Ef  field  components  (each  of  which  varies  as  JQ(Tr)with  the  same  T),  does  not  smoothly 
reduce  to  the  radial  TEM  mode  wave.  When  |t|  r  «  1,  then  JQ(Tr)  »  1  and  the  usual  one¬ 
dimensional  results  are  retrieved.  The  current  density  is  similarly  described  as 


J  +  J  =  y  J  (T.r)  cr  B.  e 
n  p  Li  o'  l  l  l 

i=l 


-jKiX 


(22) 


Fig.  2.  The  approximate  and  exact  silicon  IMPATT  diode  impedance  vs  frequency  (f) 
for  a  diode  in  which  the  avalanche  zone  (E)  is  1.5  microns  in  length,  the  drift  zone  (1)) 
is  2.75  microns,  JQ  -  lO^amp/cm2,  <yQ  -  6.6  x  103,  -  0.164.  (a)  Heal  part;  (b)  im¬ 

aginary  part. 
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REAL  part  of  z  (ohms)  REAL  PART  OF  Z  (ohms) 


FREQUENCY  (GHz)  FREQUENCE  {GHi) 

(o)  {tO 

Fig.  3.  The  approximate  and  exact  silicon  1MPATT  diode  impedance  vs  frequency  (f) 
for  a  diode  in  which  the  avalanche  zone  (L)  is  0.5  micron  in  length,  the  drift  zone  (D) 
is  0.92  micron,  J0  3.76  X  10^  amp/cm^,  cv0  2  x  104,  -  0.315.  (a)  Heal  part, 

(b)  imaginary  part. 


(a)  (b) 

P  ig.  4.  The  approximate  and  exact  silicon  I MPATT  diode  impedance  vs  frequency  (f) 
for  a  diode  in  which  the  avalanche  zone  (L)  is  0.32  micron  in  length,  the  drift  zone  (1)) 
is  0.53  micron,  J0  7.18  x  lO^amp/cm^,  cvD  3.125  x  10^,  Oq  0.427.  (a)  Heal  par-1, 

(b)  imaginary  part. 
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REAL  PART  OF  Z  (ohms) 


(o) 


50  60  TO 


FREOUENCY  (GHr) 
(b) 


Fig.  5.  The  approximate  and  exact  silicon  IMPATT  diode  impedance  vs  frequency  (f) 
for  a  diode  in  which  the  avalanche  zone  (L)  is  0.22  micron  in  length,  the  drift  zone  (D) 
is  0.387  micron,  JQ  =  1.035  X  104  amp/cm2,  rv0  =  4.545  X  104,  a'Q  -  0.542.  (a)  Real 

part;  (b)  imaginary  part. 


Fig.  6.  The  real  and  imaginary  parts  of  T^  vs  frequency  (f)  for  an  X-band  diode  with 


1.5  microns,  JQ  =  10^  amp/cm2,  oQ  6.6 


X  10- 


-  0.164,  and  *  -  10"^  farads/cm. 


U  - r  ' -  '  U  '  U  '  T 

The  square  of  the  radial  wavenumber  (for  a  TEM  cylindrical  wave)  K#  (which  is  purely 
real)  is  indicated  for  comparison,  (a)  Re  {T2};  (b)  Im  {T2}. 
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APPENDIX 

IMPEDANCE  DERIVATION 


'Phe  continuity  equations  for  electron  and  hole  currents  are^ 
On/cM  ( l/q)  a.ln/dx  +  g  , 

Op/at  —  (i /q )  a.ip/dx  +  g  , 


(23) 

(24) 


where  g  nQ\'s(n  f  p)  4  H^(nQ  4  Pt))  is  the  small-signal  AC’  component  of  the  electron-hole 

generation  rate,  Y  is  the  saturated  drift  speed  of  the  carriers,  n  and  p  are  the  AC  electron 

and  hole  currenl  densities,  and  n  +  p  are  the  DC  electron  and  hole  currenl  densities.  After 

o  o 

some  algebraic  manipulation,  along  with  the  assumption  that  all  quantities  vary  as  e,,a"  ,  Kqs.{23) 
and  (24)  may  be  written 


1 )  n  —  o-  V  p  -  C  E  0 
n  o  s  ^  x 


D  p  -  (v  V  n  —  C  E  0 
p^  o  s  \ 


(2S) 

(2M 


whe  re 


D  _  (a/ot)  -  Y  r)  /f) x  —  a  V  ) 
n  s'  os 


Dp  _  (f)/9t  4  Vg<>/0x  —  Of q \  g ) 


C  o'  \  ^(n  4  p  ) 

o  s  o  o 


From  these  equations  it  can  be  readily  shown  that 


[D  1)  (o'  V  rj  J  -  qV  C  [D  +  a  V  ]  E 

lnp  o  s  1  n  ns  1  p  os'x 


[D  D  -  (or  V  )  1  J  qV  C  [D  4  o  \  ]  h 

1  n  p  o  s  1  p  n  s  1  n  os'  x 


(^7 ) 


jU’t-K.x) 

For  e  dependence  it  is  recognized  that  Eqs.  (6),  (23),  and  (24)  yield  a  third-order  sys¬ 

tem  of  equations  for  which  the  euhic  dispersion  relation  reveals  the  three  roots  for  K.  given  in 

Sec.  II.  Then  E  ,  .1  ,  and  .1  nnisl  have  the  form 
x  n  p 


V 


E. 


•tKix  \ 
n.  e  ) 


(71 


3 

V 

i  1 


a  .  E. 

ni  l 


V 


a  IT 


p  — <  ]^ii 

i  1 


(2K) 


where  ct  .  and  a  .  can  be  determined  hv  substituting  Kq.  (7)  into  Eqs.  (27)  and  turn  out  to  be  given 
m  pt 

by 


CTi  1 

r  VKi 

(J  .  — 

i  +  — - 

ni  2 

1 

r  VK 

(T  .  — 

1  - - 

Pi  2 

l 

( 29 ) 


i  Eq. 

(9). 

Applying  the  usual  boundary  conditions 

-  0 

at 

x  =  -L/2 

-  0 

at 

x  L/2 

(30) 

it  follows  that  B^  B^  and 


b2/b1  =  B3/Bl 


M 


(31) 


where  M  is  defined  in  Eq.  (17).  Thus,  Eq.(7)  can  be  written 
Ex  -  B^j  [1  +  M  cos  (K  X)] 

Eq.  (18)  is  obtained  by  summing  Eqs.  (28). 
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